SIMULTANEOUS ESTIMATES FOR VECTOR- VALUED GABOR 
FRAMES OF HERMITE FUNCTIONS 

H. FUHR 

Abstract. We derive frame estimates for vector-valued Gabor systems with window 
functions belonging to Schwartz space. The main result provides frame bound estimates 
for windows composed of Hermite functions. The proof is based on a recently established 
sampling theorem for the simply connected Heisenberg group, which is translated to a 
family of frame estimates via a direct integral decomposition. 



1. Introduction and main results 

The aim of this paper is to derive frame estimates for vector-valued Gabor systems. We 
consider the space L 2 (M; C d ) of vector- valued signals. Elements of this space can also be 
understood as vectors f = (fi, ■ ■ ■ , fd), with f t G L 2 (IR; C), which amounts to identifying 
L 2 (M; C d ) with the d-fold direct sum of L 2 (M). Gabor systems in this space are obtained 
by picking a window function f G L 2 (IR; C d ), and applying translations and modulations. 
Here the translation and modulation operators are given as 

(T y f)(x) = f(x - y) , (M c f)(x) = e 2 ^f(x) (y,£ e R) 

For any 7 = (71,72) G M 2 , we denote the associated time-frequency shift of a function 



f G L 2 (M; C d ) by 



Now, given a lattice T C M 2 and f G L 2 (M; C d ), the resulting Gabor system Q(f,T) is 
given by 

£(f,r) = (f 7 ) 76r . 

Recall that a family (f]i)iei of vectors in a Hilbert space TC is called a frame if it satisfies 

(i) am 2 <J2\M\ 2 < £|MI 2 , 

16/ 

for all <£> 6 W, with constants < A < 5. These constants are called frame bounds; they 
are generally not assumed to be optimal. A Gabor system Q(f, T) that is also a frame is 
called a Gabor frame. 

For the case c? = 1, Gabor frames have been studied extensively; see e.g. ^3]- For a 
treatment of the case d > 1, in somewhat different terminology, see e.g. 

The problem of constructing a Gabor frame in dimension d contains that of constructing 
d Gabor frames: If Q(f, T) is a frame with bounds A, B, then Q(fi,T) is a frame with 
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frame bounds A, B, for each component f\ of f. More generally, Q{f,T) is a frame of 
L 2 (IR; C d ), for f = (f\, . . . , fj), with d < d, again with frame bounds A and B. 

But the converse need not be true: Note that by definition of the scalar product in 
L 2 (M;C d ), 

d 

(g,fx) = ^2(9i,fi\x) , 

i=l 

where we used fi\\ to denote the action of a time-frequency shift on fi. Hence, can- 
cellation may prevent the higher-dimensional system from being a frame even when all 
components fi, ■ ■ ■ , fd generate a frame. In fact, one can easily see that a obvious neces- 
sary requirement for f to generate a frame is linear independence of its entries fi, ■ ■ ■ , fd- 
Observe however, that at least the upper frame bound for f can be estimated in terms 
of the upper frame bounds for the f: If B , B 1 , . . . , B& are optimal upper framebounds 
for f, fi, . . . , fd, respectively, then the Cauchy-Schwartz inequality entails 



(2) B<dJ2B, 



i=l 



Probably the most-studied window function for the case d = 1 has been the Gaussian, 
g(x) = Tx^ l / A e~ x I 2 . This is partly due to historical reasons: Gabor suggested using Gauss- 
ian windows JU], and the characterization of densities for Gabor frames with Gaussian 
window took more than 30 years to be fully clarified flj El • The choice of this window 
function is motivated by the way Gabor systems are employed: They are designed to 
measure time-frequency content in a signal. By the Heisenberg uncertainty relation a 
Gaussian window has optimal time-frequency concentration, and thus can be expected 
to yield a good time-frequency resolution. Moreover, for the Gaussian window powerful 
tools from complex analysis can be employed to study sampling [HJ HH] , which adds to 
its theoretical appeal. 

In this paper, we intend to derive frame estimates for window functions f consisting of 
the first d + 1 Hermite functions. For n G N , we define the n th Hermite function h n by 



— I)" „.2 in d 



2 



(3) hJx) = ; ^ e x 12 ^e~ x . 

where the normalization factor ensures ||^ n ||2 = 1- The above defined Gaussian equals 
Hq, whence the problem considered here can be viewed as a generalization of Gabor's 
original question. 

The vector- valued windows that we are interested in are given by 

h d = Cio,.,MeL 2 (R;C w ). 

We intend to give frame conditions and estimates for Q(h d , Ai(Z 2 )), with Ai an invertible 
matrix, in terms of a matrix norm defined by 

(4) ||M|| =sup{||^|| 2 : IML < 1/2} . 

Here ||z|| p denotes the usual £ p -norm on M 2 . This choice of matrix norm may seem 
somewhat peculiar, and in fact the theorems below can be formulated with respect to 
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any other norm on matrix space. The use of (@J emphasizes the close connection to 
sampling estimates on the Heisenberg group. 

Let us now state the main results of this paper. The chief purpose of Theorem 11.11 is 
to allow a better understanding and formulation of Theorem 11.21 Nonetheless, Theorem 
11.11 is of independent interest. Even though we expect it somehow to be part of Gabor 
analysis folklore, we are not aware of any previous source for this result; not even for 
d=l. 

Theorem 1.1. Let f = (f 1} . . . , f d ) e L 2 (R, C d ) be given, with fa e S{R), and (f i} fj) = 
f or oil 1 < i,j < d. Then there exists a constant < Cf < 1 such that for all 
matrices M. with ||.M|| < Cf, the system Q(f, M.{1?)) is a frame with frame constants 

|det(A)| \ L + C t ) ■ 

Hence the tightness of the frame estimate, which is the quotient of the two frame 
bounds, approaches 1 as — > 0, with speed proportional to \\M. \\ 2 /C^. 

We observe that Theorem 11.11 holds for the supremum Cf of all possible constants, and 
that this choice provides the sharpest possible statement. Then the main result of this 
paper is the following estimate: 

Theorem 1.2. There exists a constant Cm < 1 such that for all d e No 
(5) C* hd > °^ . 

2. Method of proof: Sampling vs. frame estimates 

The proof of Theorem 11.21 relies on a combination of various techniques: It uses a 
sampling estimate for the Paley- Wiener space PW(M), established in jU]. The definition 
of this space uses a particular differential operator on HI, the so-called sub-Laplacian. 
Hermite functions enter in the spectral decomposition of this operator, and it is this 
connection that will allow to relate the sampling theorem to frame estimates for Hermite 
functions. 

The connection between frames and sampling theory is not exactly new, in fact it is 
at the base of frame theory, which originated from nonharmonic Fourier series and their 
connections to irregular sampling over the reals, see the classic paper |3] . For the sake of 
explicitness, assume we are given a sequence A = (\k)k& °f sampling points in M. We 
are looking for criteria that allow to reconstruct Paley- Wiener functions, i.e. / G L 2 (IR) 
whose Fourier transform has support in the unit interval [—0.5,0.5], in a stable manner 
from its restriction to A. Noting that 

/(A*)= f " m)e 2niXk ^=(f,e Xk ) , 

J -0.5 

we find the following two equivalent conditions, with identical constants A and B in both 
cases: 

(1) The sequence A fulfills the sampling estimate 



A\\f\\i<Y,\f^)\ 2 < B 



2 

2 ) 
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for all Paley- Wiener functions. 
(2) The sequence (e\ k )k£Z fulfills the frame estimate 

A\\F\\l<J2me Xk )\ 2 <B\\F\\l, 

feez 

for all F G L 2 ([-0.5,0.5]). 

This equivalence can be used in two ways: For instance, observing that the choice 
Afc = k [k G Z) results in the Fourier orthonormal basis [tk)k&L of L 2 ([— 0.5, 0.5]), the 
implication (2) =^> (1) leads to Shannon's sampling theorem. Conversely, for irregular 
sampling sets, condition (1) can often be checked using tools from complex analysis tech- 
niques, and then (1) =^> (2) results in frame estimates for irregularly spaced exponentials. 

In this paper, we use a similar approach for the Heisenberg group H: This time, 
previously established sampling estimates for Paley- Wiener functions on EI will allow to 
derive frame estimates for Hermite functions, by an analogue of the implication (1) =>- (2). 
For this purpose we will need to work out the connections between PW(HI) and the 
Hermite functions, which is the topic of Sections 4 and 5. But first let us prove Theorem 

o 



3. Proof of Theorem 11.11 

Let Vf : L 2 (M; <C d ) — > L 2 (M 2 ) denote the windowed Fourier transform, 

V f g(x,0 = (g,T x M^) . 

The orthogonality relations for the windowed Fourier transform and the pairwise orthog- 
onality of the components of f imply that Vf is an isometry. Hence its image 7if is a 
closed subspace of L 2 (M 2 ). As outlined in the previous subsection, the isometry property 
of Vf implies that a frame estimate for Q(f, M.(J?)) is the same as a sampling estimate 
for 7Yf, with sampling set M.{1?). 

We intend to utilize the techniques from for this purpose, hence we will need oscil- 
lation estimates. We define 

0SCr(/)(ar) = sup \f(x) - f(y)\ . 

\x-y\<r 

Our first aim is to prove 

(6) ||osc r (F)|| 2 <^-||F|| 2 ,VFGWf , 

Of 

for a suitable constant Cf. 

For this purpose we first observe that the projection Pf onto Tif is obtained by twisted 
convolution [T3j: 

(P f G)(x,0 = 

= / G(x', g)F{x - x',£ - Oe^'-^Wf^' , 

where we let F = Vff . 
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Hence G = G$F for G G Hf, and therefore 
osc r (G)(x,£) 



sup 

|(x,e)-(a«"4")|<r 



I G(x', (F(x -x',t- O - F(x" - x', C - O) e^'-^dx'dt 1 
< [ \G(x',0\ sup \ F {x-x'^-0-F{x" -x',e ~0\dx'di' 

Jw? \(x£)-(x",£")\<r 

= [ \G(x',£)\osc r (F)(x-x',£-C)dx'd£ 
= \G\ *osc r (F) , 

where convolution is taken with reference to the abelian group structure on R 2 . Hence 
(7) ||osc r (G)|| 2 < UGHallosc^F)!!! . 

Now the second factor can be estimated by 



|osc r (F)||] 



sup \F(x,£)-F(x?,Z')\dxdt 



< I r sup 



dF 



+ 



OF 



dxd£ 



< r 



I V \\D a F\\^ Br{x ^dxdi 

J ^ 2 \a\<l 

[ V ||D tt F|| liBr(a , €) da;d£ 



< r 

|q|<4 

r ^ I^ICbJI^FH! 

|a|<4 



where we used the mean value theorem for the first inequality, and the Sobolev embedding 
theorem for the third. Here Cs r denotes the norm of the embedding W 3,l (B r ) <^-> C(B r ) 
PU Theorem 5.4]. Clearly, \B r \ = r 2 \Bi\; on the other hand, a dilation argument estab- 
lishes for r < 1 that Cn < t~ 2 Cb x - Hence 



osc 



.{F^KrlB^CB, \\ D * F \ 

\a\<4 



I ■ 



and thus (J7J) implies (JBj), with Cf 



\B x \C Bl J2 W DaF \ 

|a|<4 



Now by letting K = M([— 0.5, 0.5] 2 ), we obtain M 2 = U 7 ex(z 2 ) 7 + -ft' as a disjoint 
union, and K has Lebesgue measure |det(.M)|. Moreover, by definition of the norm, 
K C B r , for r = ||.M||. Now, by [HI Theorem 3.5], the oscillation estimate © results in 
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the sampling estimate 

^ i 1 - k )' ^ £ X, |F(7)|2 ~ m i 1 + i ) v »* ■ 

which is Theorem ll.il □ 

Remark 3.1. As the proof shows, one can weaken the requirement fa G <S(R) to the 
membership of Vff in a suitable Sobolev space. In addition, pairwise orthogonality of the 
/, can be weakened to linear independence; in this case the frame estimate will contain a 
term involving the Gramian matrix of the /j. □ 



4. Fourier transform on the Heisenberg group 

The (simply connected) Heisenberg group is defined as EI = R 3 , with group law 

(P, Q, t)(p, q', t') = (p + p , ,q + q\ t + t'+ (pq - p'q)/2) . 

For the following facts concerning H, we refer the reader to jHj- H is a step-two nilpotent 
Lie group, with center Z{M) = {0} x {0} x R. EI is unimodular, with two-sided invariant 
measure on H given by the usual Lebesgue measure of R 3 . 

Given A G R* = R \ {0}, the Schrddinger representation p\ of EI acts on L 2 (R) via 

p x (p,q,t)=e- 2 ^ t - p ^T Xp M q . 

This is an irreducible unitary representation of EI. The family of Schrddinger represen- 
tations provides the basis for the Plancherel transform of the group, a tool that is of key 
importance for this paper. 

Before we describe this transform in more detail, let us quickly recall the basics of 
Hilbert-Schmidt operators: The space of Hilbert-Schmidt operators on a Hilbert space H 
is given by all bounded linear operators such that 

\\n 2 HS = y: \\Tvi\\ 2 

iel 

is finite; here (?7i)ie/ denotes an arbitrary orthonormal basis of 7i. It is well-known that 
the norm is independent of the choice of basis, and defines a Hilbert space structure on 
the Hilbert-Schmidt operators with scalar product 

(S,T) = trace(T*S) = ^<%,T^> . 

Now, given / G L 1 H L 2 (H), one defines 

Pxif) = / P\(x)f(x)dx , 
Jw 

understood in the weak operator sense. This is just the canonical extension of the rep- 
resentation pa to the convolution algebra L 1 (H). The mapping / i— > (px(f))xm* is called 
the group Fourier transform. This nomenclature is justified by the observation that the 
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euclidian Fourier transform is obtained by integration against the characters of the addi- 
tive group. Moreover, it turns out that for / G L 1 PI L 2 (H), the group Fourier transform 
is in fact a family of Hilbert-Schmidt operators, and we have the Parseval relation 

11/111= / \\Px(f)\\Hs\MdX ■ 

JR 

The Parseval relation allows to extend the Fourier transform to L 2 (H), yielding the 
Plancherel transform, a unitary map 

L 2 (H) -> / HS{L 2 {R)) \X\d\ , 
Jr 

where the right hand side denotes the direct integral of Hilbert-Schmidt spaces. We 

denote the Plancherel transform of / G L 2 (H) as / = (/(A) J . This map will play 

V / AeM* 

the same role as the euclidean Fourier transform in the discussion of Section 2. 

Of crucial importance for this usage of the Plancherel transform are its algebraic 
properties, providing a decomposition of various operators and representations acting 
on L 2 (H). For instance, if we denote the left regular representation by L x f(y) = f(x~ 1 y), 
for / G L 2 (HI) and x, y G H, then 

Px(L x f) = P\(x)o Px (f) , 

which provides the decomposition of the left regular representation L into a direct integral, 

L ~ / p A ® 1 |A|dA . 
Jr* 

Similarly, the right regular representation R decomposes by the formula p\(R x f) = Px(f)° 
Px(x)*- 

The decompositions extend to commuting operators: For any bounded operator T 
commuting with L, there exists a measurable field (Ta)a£M* of bounded operators on 
L 2 (IR) satisfying p\{T f) = p\(f) ° T\, or in direct integral notation 

T~ I Idof A \X\dX . 
Jr* 



5. Paley- Wiener space on the Heisenberg group 

In this section we outline the definition of Paley- Wiener space on H and its relation 
to Hermite functions. The central role of Hermite function in the decomposition of the 
sub-Laplacian has been observed previously, e.g. in [TT]; the results presented below can 
be found also in [T7j . 

We define a left-invariant differential operator P on EI by 

(8) (P /)(„, 9 , t ) = lim f((P^)(h,0 0))-f( M ,t) 

h— >o h 

corresponding to the subgroup R x {0} x {0}, and Q is a left-invariant operator associated 
to {0} xix {0} in the same manner. P, Q are viewed as elements of the Lie algebra f) of 
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H; we have [P, Q] = T, the infinitesimal generator of the group center. This observation 
exhibits i) as a stratified Lie algebra, 

f> = Vx + V 2 

with = span(P, Q), and V 2 = R ■ T = [V u V x ). 

Of particular interest for analysis on these groups is the sub-Laplacian; as the name 
suggests, it can be viewed as a replacement for the Laplacian over R n . For the Heisenberg 
group, this operator is defined by 

C = -P 2 - — 
4vr 2 • 

The normalization of Q is chosen for the sake of convenience. £ is a left-invariant positive 
unbounded operator on L 2 (H). We denote its spectral measure by Tic- The Paley- Wiener 
space on EI is then defined as 

P^(H)=n £ ([0,l])(L 2 (H)) . 

Note that, up to normalization, this definition is completely analogous to the definition of 
bandlimited functions on R, since the euclidean Fourier transform can also be read as the 
spectral decomposition of the Laplacian. The projection TE^ ( [0, 1]) is left-invariant, and is 
therefore decomposed by the group Fourier transform into a direct integral. The following 
lemma provides an explicit calculation of this decomposition via Hermite functions. In 
the following, we use the notation D a f(x) = \a\~ 1 / 2 f(\a\~ 1 x). As with translation and 
modulation operators, we use the same symbol for operators acting on scalar- and on 
vector-valued functions. 

Lemma 5.1. (a) (/i n )neN an orthonormal basis o/L 2 (R). 

(b) The system (/i n )neN ^ s an eigenbasis of the Hermite operator Hf(x) = x 2 f(x) — 
fix), 

(9) Hh n = (2n + l)h n . 

(c) For every real a ^ 0, the dilated system (/i n ,a)neN ; defined by 

h n>a (x) = (D\ a \i/ 2 h n )(x) = |a|~ 1/4 /i n (|a|~ 1/2 x) , 
is an eigenbasis of the scaled Hermite operator H a f(x) = x 2 f(x) — a 2 f"(x), 

(10) H a h n , a = \a\(2n + l)hn, a 

(c) The sub-Laplacian decomposes into a direct integral of scaled Hermite operators: 

p x (Cf) = Px (f) o H x , 

for all f G Q°(H). 

Proof. For part (a) confer Corollary 6.2, Theorem 6.14]. Part (b) follows from this 
by straightforward computation. Part (c) is established by formal calculation from (jHJ) and 
the analogous formula for Q, using the decomposition of the right regular representation. 
□ 
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Parts (b) and (c) contain the ingredients of the direct integral decomposition of PW(M). 
For the precise formulation of this result and its proof, the tensor product notation for 
Hilbert-Schmidt operators will be useful. Given vectors 77, ip in a Hilbert space 7i, we let 

77 <g> ip : z h-> (z, (p)rj , 

which is a rank-one operator on TC. Note that the notation is only conjugate linear in tp. 
The Hilbert-Schmidt scalar product of two elementary tensors is 

(77 ® If, 7]' ® if') HS = (7/, 7]'} n (ip', ip) n . 

Moreover, for any pair S, T of bounded operators, S o (77 <g> ip) o T = (Si]) ® (T*ip). 

Now, given any orthonormal basis (fiji^i of 7^, every Hilbert-Schmidt operator T has 
a unique decomposition 

T = ^ Vi ® • 

iei 

Hence, if S = ^ (§) </?j is another Hilbert-Schmidt operator, we obtain for the scalar 
product 



ieJ 

Observe in the formulation of the following proposition that P\ involves the first d(X) + l 
Hermite functions. For |A| > 1, we have d(X) = —1, and thus P\ = 0. 

Proposition 5.2. Letting 

d(X) = 4- - - 
L 2 |A| 2 

and 

n=0,...,d(A) 

t/ie projection onto Paley- Wiener space is given by 

(12) (H £ ([0, 1])(/)) A (A) = /(A) o P A , V/ G L 2 
T/^e operator field (P\)AeM* ^ e Plancherel transform of a function p G L 2 (H), whence 

n £ ([o J i])(/) = /*p. 

Proof. We apply the above considerations to the case = L 2 (IR) and its orthonormal 
basis (/vaWno- Hence each Hilbert-Schmidt operator T on L 2 (IR) has a decomposition 

(13) T = ^ Vn ® ^n,A , 

neN 

and we obtain from ()1U|) that 

This shows that the map T h t)„8 /i n ,A can be understood as a projection onto an 
eigenspace of the operator T 1— > T o iJ A , with associated eigenvalue |A|(2n + 1). By 
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definition of Paley- Wiener space, only eigenvalues < 1 are admitted, which shows that 
the definition of Pa indeed yields ()12|) . 

For the second statement, we compute the norm of the operator field in the direct 
integral space. First observe that P\ = for |A| > 1. Moreover, the squared Hilbert- 
Schmidt norm of a projection equals its rank, whence ||-Pa||//s = d(\) < ^kr, and thus 

/ \\P X \\ 2 HS \X\d\< f J- \x\d\=l. 

J M.* J — 1 I I 

Hence (P\)\ has a preimage p under the Plancherel transform. Finally, (|12|) and the 
convolution theorem [HI Theorem 4.18] provide that Pf — f *p. □ 
The motivation for considering PW / (H) is the existence of sampling estimates. The 
formulation of the sampling theorem requires some additional notation. We fix a quasi- 
norm | • | : H — > by 

\{p,q,t)\ = {p 2 + q 2 + \t\) l '\ 

and write B r for the unit ball around 0. A discrete subset T C H is called a quasi-lattice 
if there exists a relatively compact set K C H such that H = {J„ /f - r / yK, as a disjoint 
union. Such a set K is called complement of V. 

Theorem 5.3. [9, Theorem 5.11] There exists a constant < Ch < 1 with the following 
property: For all quasi-lattices V possessing a complement K contained in a ball of radius 
r<C m and all f e PW(U) 

(14) r^|(i - r/c H ) 2 n/i^ < i/(7)i 2 < ,4(1 + r/c m nni . 

I I 7e r I I 



6. Proof of Theorem 11.21 

We will now derive Theorem 11.21 from Theorem 15.31 basically by explicit calculation. 
The following lemma can be seen as an analog of (1) <^ (2) from Section 2. A version of 
this result was obtained in [HJ Proposition 6.11]. 

Lemma 6.1. Suppose that V C H is of the form V = F x aZ, with rcK 2 and a > 0. 
Consider the following statements: 

(a) For all f G PW{W), 



(is) A\\f\\i<j2\m\ 2 <B 



(b) For all d £ No, and for almost all A with |A| < ^TT' the system G(h d , | A 1 1//2 T) is 
a frame o/L 2 (R;C d+1 ) with frame bounds a\\\~ x A anda\\\~ 1 B. 

Then (a) =^ (b), and if a < 1/2, (b) =^ (a). 

Moreover, if T = M.{1?\ for a suitable invertible matrix, the frame estimates in (6) 
are valid for all A < • 
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(16) 



Proof. For the proof of (a) =>• (6), let / G PW / (H) be given. Then we can write 

d(A) 

/(A) = y^y»,A <8> 



%A 



i=0 



for suitable functions (p^x € L 2 (M). In the following, we also use the notations 



$A = (<A),A, ' ' 



and 



hA = (ho,x, • • • , hd(\),\) • 
Let £ C 1* be a Borel set contained in an interval / of length I /a, and consider g with 
g — f ■ lfi. Observing that g G PW / (H), we can compute the £ 2 -norm of its restriction to 
T' as follows: 

Ei^)i 2 = Ei^^mi 2 

7er' 

2 



7er' 



7er' 



(f(X),Px(j)m)\MdX 



= E 

(z,fc)er,ne 

Applying the Parseval formula for the interval, we thus obtain 



(/(A),p A (Z,M)p( 7 )>e 27nA HA|dA 



(17) 



5>WI 2 = a ~ l E / K/(A),PA(/,^0)p( 7 ))r|A| 2 rfA 
7er' {l,k)er jE 

= a- 1 I E |($A,TA*M fc |A| 1/2 h A )e" AZfc | 2 |> 
^ B («,fc)er 

= a" 1 / E |($A,T|A|iM fc |A| 1/2 h A )| 2 |A|dA 



(z,fc)er 



where the last equation used (fTTj) to express the Hilbert-Schmidt scalar products as scalar 
products of vector- valued functions, as well as symmetry of V to replace A by |A|. 
On the other hand, by the Plancherel formula, we find 



I^III 



\\f{\)\\ 2 H S \m 



I^A|| L 2( KiC d(A)) 



X\d\ . 



Hence the lower sampling estimate yields 
> II 2 I 

^llL 2 (IR,C d ( A ))l 



IN) A I IH'^^^lMdX^a- 1 ! E \{§x,T w M k \\\ l / 2 h x )\ 2 \\\d\ 

Je (i,k)er 
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Since this inequality holds true for all Borel sets E of diameter at most 1/a, it has to 
hold pointwise a.e. for the integrands, i.e. after shifting constants: 

(19) a\\\- x A\\® x \\l^ cd{X)) < \(<5>x,T lMl M k h x )\ 2 (a.e.A) . 

(/,k)er 

This is already quite close to the desired lower frame estimate, except that it holds on a 
set of A's which may depend on the choice of / (or equivalently, on the field ($a)a£K*)- 

The next step is to establish (J19)) for all / G PW(EF) and all A in a set with complement 
of measure zero, independent of /. For this purpose we pick a sequence (/ n )neN with dense 
span in PW(M), and obtain a set Q C [—1, 1] with complement of measure zero such that 
(jig) holds for all A G tt and all / in the (Q + iQ)-span of (/ n ) neN . But then the (Q+iQ)- 
span of (f(X))n&i * s dense in HS(L 2 (R)) o P A , for all A in a Borel set Q' C [—1, 1] with 
complement of measure zero. Hence for all A G Q fl Q', the frame estimate holds on a 
dense subset of HS(L 2 {R)) o P A , which is sufficient. 

Thus we have finally established (JED for almost all A G [-1,1], and all $ A G L 2 (M, C d(A) ). 
The same argument applies to show the upper estimate with constant a|A| _1 P. 

Now, using h x = P| A |i/ 2 h d ( A ), and the relations 

M^D b = D b M bi , T x D b = D b T b -i x 

we find that 

T w M k h x = D w/2 (T w/2l M w/2k h d ^) . 

Since the image of a frame under a unitary map is a frame with identical constants, we 
finally obtain that (T^ x ^/2iM^ x ^/2 k h d ^)(i : f i .- )& Y is a frame, for almost all |A| < 1. Now part 
(b) follows from d < d(X) for A < 

For the converse direction observe that by assumption on a, all Plancherel transforms 
of elements of PW(HI) are supported in an interval of length 1/a. Hence (JT7j) holds for 
all g G Piy(IHf), where this time E = [—1,1], and the field ($a)a corresponds to the 
Plancherel transform of g. But then (b) =^ (a) is immediate. 

The proof that the "almost everywhere" contained in the statement can be omitted for 
lattices relies on semi-continuity properties of the frame bounds. 

For any unit vector f G L 2 (IR, C d ), consider the function 

© f :]0,^ T [3A^ l(f^A4|Arh d >| 2 . 
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We compute 

©f(A) = |A| l(f m ,^M fc h d )| 2 

(l,k)e\\\V 2 M(Z 2 ) 

d 

(/,fc)e|A| 1 / 2 A4(s 2 ) i=o 
= l A l EE<^' TiMkh^TiMkhi) 

i,j (l,k) 

= \x\J2<^J2^ TiMkh ^ TiMkhi ) 

id (l,k) 

= \M ^2{fh S haii .\ X \i/2 M fj) . 

id 

Here we used the linear operator S gitgr ^ associated to functions gi,<?2 and an invertible 
matrix Af, defined by 

SgiM-Af) = E TiM k gi)TiM k g2 ■ 

(i,k)eAf(z 2 ) 

Theorem 3.6] states that S : M 1 (R) x M 1 (E) x GL(2,R) -> £(L 2 (R)) is continu- 
ous, where the right-hand side denotes the space of bounded operators endowed with 
the norm topology, and M 1 (M.) is the Feichtinger algebra; see e.g. for a defini- 

tion and basic properties. Now the inclusion <S(R) C M 1 (R) entails that the map 
A i— > (/j, S huhj:> \ X \i/2Mfj) is continuous, for all < i, j < d, and then f is continuous. 
Next consider the map associating to each A the optimal upper frame bound, given by 

Bopt :]0, ttt—^ A i-> sup f (A) . 
ztz + i ||f ll=i 

The supremum is always finite: By |T3*1 Corollary 6.2.3], the frame operator of a one- 
dimensional window in the Schwartz class is always bounded. Hence the upper frame 
bound also exists in the vector valued case, by (J2J). 

As the supremum of a family of continuous functions, B opt is lower semi-continuous, and 
then A i— > _B op t(A)|A|~ 1//2 is lower semi-continuous as well. We already know that the latter 
map is bounded from above by aB on subset of ]0, with complement of measure 
zero. This subset is dense, hence lower semi-continuity implies -B op t(A)|A| -1 / 2 < aB on 
the whole interval. 

The analogous reasoning, replacing lower by upper semi-continuity, applies to the lower 
frame bound, and we are done. □ 

Proof of Theorem 11.21 Fix d G N . Suppose that Ai is given with \\Ai\\ < 
C m /V2d+1. Let K = M([-0.5, 0.5) 2 ), then if is a complement of M(Z 2 ) in R 2 , con- 
tained in a ball of radius tq = ||-M||, and with measure |det(A / f)|. Moreover, by choosing 
a > small enough, the set K' = K x [—a/2, a/2) is contained in a ball of radius 
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r + e < C^/\/2d + 1 (with respect to the quasi-norm on H). In addition, K' is a com- 
plement of r" = M.(l?) x aL: For any (p,q,t) G H, there exist unique (pi,gi) G T and 
(P2, 92) G if with (px +P2, 9i + 92) — {Pi q), an d finally unique I G Z and s G [—a/2, a/2) 
with s + al = t — {pxq% — p2?i)/2. But these choices imply (pi, g 1; a/)(p 2 , g 2 , s) = (p, g, t). 

We will apply Theorem 15.31 to dilated copies of T'. For this purpose, let for a > 
and (p,q,t) G HI, 5 a (p,q,t) = (ap,aq,a 2 t). It is easy to check that this defines a group 
automorphism 5 a fulfilling \5 a (p, q, t)\ = a\(p,q,t)\. Hence S a (T') is a quasi-lattice, with 
complement S a (K) contained in a ball of radius a(r + e). Hence, for any a < Ce/(r + e), 
the sampling theorem provides the estimate 

1 A a(r + e) x " 



a 4 |det(A^)|a 



765 a (r') 

< 1 ft , a(^o + ey 



2 

2 ■ 



a 4 |det(.M)|a \ Ce 
An application of Lemma f6. II then yields, for all A < 2^r> that ^(h d , alAI 1 / 2 ^ (Z 2 )) is a 
frame with bounds ffl 2| A | |d c t(A^)| {} "F °^ Cb!'^ ) ■ Letting a 2 |A| = 1 provides a lower frame 



bound for Q(h d , M. (Z 2 )) given by 
sup 



{i5sk( i - !i ^) ,: ^ fTT ^<^ + 4- 

Observe that the restriction a > \/2d + 1 is imposed by |A| < • By monotonicity, the 

supremum is ^ dct ^ M ^ fl — v/2d+ ^ r ° + — j . Sending e to zero provides the lower estimate of 
Theorem 11.21 The upper estimate is obtained in the same fashion. □ 



Concluding remarks 

It is a standard observation that the construction of Gabor frames is equivalent to the 
discretization of the inversion formula associated to a certain discrete-series representation 
of the so-called reduced Heisenberg group M r , which is the quotient of HI by a discrete 
central subgroup; see e.g. [TH] . The proof of Theorem 11.21 shows that in working with 
H one needs to deal with a fair amount of additional technical details (in particular due 
to the occurrence of direct integrals), that one avoids by considering Et r . The benefit of 
this approach lies in the fact that a single sampling estimate, namely ()14j) . gives rise to 
a whole family of Gabor frame estimates, namely (|5]L valid for all d > 0. 

The main results of this paper provide rather intuitive asymptotic estimates for Ga- 
bor frame bounds. A major drawback of these estimates is that they involve unknown 
constants. A "formula" for C*e is given in [3], involving operator norms for differential 
operators on Fl^(H) as well as a Sobolev constant for the unit ball H; the argument 
is very similar to the estimate of the constant Cf in the proof of Theorem 11.11 While 
rough estimates for the differential operators should be obtainable from the Plancherel 
transform, which decomposes the differential operators as well as PW(M), we are not 
aware of a reasonable estimate for the Sobolev constant for HI. In any case, we stress 
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that the constant Ch in the sampling theorem is the same as in Theorem 11.21 this was 
the chief motivation for picking the matrix norm (J3J). 

For single Hermite functions, the results obtained here compare in an interesting way 
with recent results due to Grochenig and Lyubarskii. Using complex analysis methods, 
they obtained the following statement [T2"l Theorem 3.1]: 

Theorem 6.2. If | det M | < (d + thenQ{h d ,M) is a frame for L 2 (R) . 

For the isotropic case, i.e., M. = a ■ Id, this result provides a criterion that is very close 
to our Theorem 11.21 Any a below a threshold ~ rr 1 ! 2 guarantees a frame. In the general 
case however, Theorem 16.21 is much more widely applicable: At the same time | det A4 \ 
can be made arbitrarily small and ||.M|| arbitrarily large. 

On the other hand, Theorem 16.21 does not provide frame bound estimates, and it only 
applies to the scalar-valued case. 

Let us finally comment on possible generalizations. The first possible extension consists 
in replacing R by M. n , i.e. studying vector-valued Gabor frames in L 2 (M ra ; C d ). One now 
considers the 2n + 1-dimensional Heisenberg group H n . This is a stratified Lie group, 
possessing a sub-Laplacian, Paley- Wiener space and, finally, a sampling theorem jHj. As 
for the one- dimensional case, the spectral decomposition of the sub-Laplacian involves 
Hermite functions, and an adaptation of the arguments for H should be a straightforward 
task, somewhat aggravated by additional bookkeeping. 

A second, more interesting but also more challenging type of generalization concerns the 
sampling sets, which could also be irregular. There already exists an irregular sampling 
theorem for H, however, in the transfer of the associated sampling estimates to Gabor 
frame estimates, we are crucially relying on the lattice structure of the sampling set. In 
this context, the key result is the continuity statement jSJ Theorem 3.6], and the proof of 
this result makes full use of Gabor theory developed for lattices. 

As a result, we can currently only prove statements of the following form: For all d G No 
and all uniformly discrete and uniformly dense sets T C M 2 there exists a range (0, a^) of 
dilation parameters such that Q(h d ,aT) is a frame, for almost all a 6 (0,6^), including 
an estimate of the frame bounds. Moreover, the threshold is of the order d^ 1 ^ 2 . 
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